Scalable neighbour search and alignment with uvaia.

Despite millions of SARS-CoV-2 genomes being sequenced and shared globally, manipulating such data sets is still challenging, especially selecting sequences for focused phylogenetic analysis. We present a novel method, uvaia, which is based on partial and exact sequence similarity for quickly extracting database sequences similar to query sequences of interest. Many SARS-CoV-2 phylogenetic analyses rely on very low numbers of ambiguous sites as a measure of quality since ambiguous sites do not contribute to single nucleotide polymorphism (SNP) differences. Uvaia overcomes this limitation by using measures of sequence similarity which consider partially ambiguous sites, allowing for more ambiguous sequences to be included in the analysis if needed. Such fine-grained definition of similarity allows not only for better phylogenetic analyses, but could also lead to improved classification and biogeographical inferences. Uvaia works natively with compressed files, can use multiple cores and efficiently utilises memory, being able to analyse large data sets on a standard desktop.

Repurposing an integrated national influenza platform for genomic surveillance of SARS-CoV-2 in Ghana: a molecular epidemiological analysis.

BACKGROUND: Genomic surveillance of SARS-CoV-2 is crucial for monitoring the spread of COVID-19 and guiding public health decisions, but the capacity for SARS-CoV-2 testing and sequencing in Africa is low. We integrated SARS-CoV-2 surveillance into an existing influenza surveillance network with the aim of providing insights into SARS-CoV-2 transmission and genomics in Ghana. METHODS: In this molecular epidemiological analysis, which is part of a wider multifaceted prospective observational study, we collected national SARS-CoV-2 test data from 35 sites across 16 regions in Ghana from Sept 1, 2020, to Nov 30, 2021, via the Ghanaian integrated influenza and SARS-CoV-2 surveillance network. SARS-CoV-2-positive samples collected through this integrated national influenza surveillance network and from international travellers arriving in Accra were sequenced with Oxford Nanopore Technology sequencing and the ARTIC tiled amplicon method. The sequence lineages were typed with pangolin and the phylogenetic analysis was conducted with IQ-Tree2 and TreeTime. FINDINGS: During the study period, 5495 samples were submitted for diagnostic testing through the national influenza surveillance network (2121 [46·1%] of 4021 samples with complete demographic data were from female individuals and 2479 [53·9%] of 4021 samples were from male individuals). We also obtained 2289 samples from travellers who arrived in Accra and had a positive lateral flow test, of whom 1626 (71·0%, 95% CI 69·1-72·9) were confirmed to be SARS-CoV-2 positive. Co-circulation of influenza and SARS-CoV-2 in Ghana was detected, with increased cases of influenza in November, 2020, November, 2021, and January and June, 2021. In 4124 samples from individuals with influenza-like illness, SARS-CoV-2 was identified in 583 (14·1%, 95% CI 13·1-15·2) samples and influenza in 356 (8·6%, 7·8-9·5). Conversely, in 476 samples from individuals with of severe acute respiratory illness, SARS-CoV-2 was detected in 58 (12·2%, 9·5-15·5) samples and influenza in 95 (19·9%, 16·5-23·9). We detected four waves of SARS-CoV-2 infections in Ghana; each wave was driven by a different variant: B.1 and B.1.1 were the most prevalent lineages in wave 1, alpha (B.1.1.7) was responsible for wave 2, delta (B.1.617.2) and its sublineages (closely related to delta genomes from India) were responsible for wave 3, and omicron variants were responsible for wave 4. We detected omicron variants among 47 (32%) of 145 samples from travellers during the start of the omicron spread in Ghana (wave 4). INTERPRETATION: This study shows the value of repurposing existing influenza surveillance platforms to monitor SARS-CoV-2. Influenza continued to circulate in Ghana in 2020 and 2021, and remained a major cause of severe acute respiratory illness. We detected importations of SARS-CoV-2 variants into Ghana, including those that did or did not lead to onward community transmission. Investment in strengthening national influenza surveillance platforms in low-income and middle-income countries has potential for ongoing monitoring of SARS-CoV-2 and future pandemics. FUNDING: The EDCTP2 programme supported by the EU.

Replacement of the Alpha variant of SARS-CoV-2 by the Delta variant in Lebanon between April and June 2021.

The COVID-19 pandemic continues to expand globally, with case numbers rising in many areas of the world, including the Eastern Mediterranean Region. Lebanon experienced its largest wave of COVID-19 infections from January to April 2021. Limited genomic surveillance was undertaken, with just 26 SARS-CoV-2 genomes available for this period, nine of which were from travellers from Lebanon detected by other countries. Additional genome sequencing is thus needed to allow surveillance of variants in circulation. In total, 905 SARS-CoV-2 genomes were sequenced using the ARTIC protocol. The genomes were derived from SARS-CoV-2-positive samples, selected retrospectively from the sentinel COVID-19 surveillance network, to capture diversity of location, sampling time, sex, nationality and age. Although 16 PANGO lineages were circulating in Lebanon in January 2021, by February there were just four, with the Alpha variant accounting for 97 % of samples. In the following 2 months, all samples contained the Alpha variant. However, this had changed dramatically by June and July 2021, when all samples belonged to the Delta variant. This study documents a ten-fold increase in the number of SARS-CoV-2 genomes available from Lebanon. The Alpha variant, first detected in the UK, rapidly swept through Lebanon, causing the country's largest wave to date, which peaked in January 2021. The Alpha variant was introduced to Lebanon multiple times despite travel restrictions, but the source of these introductions remains uncertain. The Delta variant was detected in Gambia in travellers from Lebanon in mid-May, suggesting community transmission in Lebanon several weeks before this variant was detected in the country. Prospective sequencing in June/July 2021 showed that the Delta variant had completely replaced the Alpha variant in under 6 weeks.

Dynamics of competing SARS-CoV-2 variants during the Omicron epidemic in England.

The SARS-CoV-2 pandemic has been characterised by the regular emergence of genomic variants. With natural and vaccine-induced population immunity at high levels, evolutionary pressure favours variants better able to evade SARS-CoV-2 neutralising antibodies. The Omicron variant (first detected in November 2021) exhibited a high degree of immune evasion, leading to increased infection rates worldwide. However, estimates of the magnitude of this Omicron wave have often relied on routine testing data, which are prone to several biases. Using data from the REal-time Assessment of Community Transmission-1 (REACT-1) study, a series of cross-sectional surveys assessing prevalence of SARS-CoV-2 infection in England, we estimated the dynamics of England's Omicron wave (from 9 September 2021 to 1 March 2022). We estimate an initial peak in national Omicron prevalence of 6.89% (5.34%, 10.61%) during January 2022, followed by a resurgence in SARS-CoV-2 infections as the more transmissible Omicron sub-lineage, BA.2 replaced BA.1 and BA.1.1. Assuming the emergence of further distinct variants, intermittent epidemics of similar magnitudes may become the 'new normal'.

SARS-CoV-2 lineage dynamics in England from September to November 2021: high diversity of Delta sub-lineages and increased transmissibility of AY.4.2.

BACKGROUND: Since the emergence of SARS-CoV-2, evolutionary pressure has driven large increases in the transmissibility of the virus. However, with increasing levels of immunity through vaccination and natural infection the evolutionary pressure will switch towards immune escape. Genomic surveillance in regions of high immunity is crucial in detecting emerging variants that can more successfully navigate the immune landscape. METHODS: We present phylogenetic relationships and lineage dynamics within England (a country with high levels of immunity), as inferred from a random community sample of individuals who provided a self-administered throat and nose swab for rt-PCR testing as part of the REal-time Assessment of Community Transmission-1 (REACT-1) study. During round 14 (9 September-27 September 2021) and 15 (19 October-5 November 2021) lineages were determined for 1322 positive individuals, with 27.1% of those which reported their symptom status reporting no symptoms in the previous month. RESULTS: We identified 44 unique lineages, all of which were Delta or Delta sub-lineages, and found a reduction in their mutation rate over the study period. The proportion of the Delta sub-lineage AY.4.2 was increasing, with a reproduction number 15% (95% CI 8-23%) greater than the most prevalent lineage, AY.4. Further, AY.4.2 was less associated with the most predictive COVID-19 symptoms (p = 0.029) and had a reduced mutation rate (p = 0.050). Both AY.4.2 and AY.4 were found to be geographically clustered in September but this was no longer the case by late October/early November, with only the lineage AY.6 exhibiting clustering towards the South of England. CONCLUSIONS: As SARS-CoV-2 moves towards endemicity and new variants emerge, genomic data obtained from random community samples can augment routine surveillance data without the potential biases introduced due to higher sampling rates of symptomatic individuals.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages.

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

Datasets for benchmarking antimicrobial resistance genes in bacterial metagenomic and whole genome sequencing.

Whole genome sequencing (WGS) is a key tool in identifying and characterising disease-associated bacteria across clinical, agricultural, and environmental contexts. One increasingly common use of genomic and metagenomic sequencing is in identifying the type and range of antimicrobial resistance (AMR) genes present in bacterial isolates in order to make predictions regarding their AMR phenotype. However, there are a large number of alternative bioinformatics software and pipelines available, which can lead to dissimilar results. It is, therefore, vital that researchers carefully evaluate their genomic and metagenomic AMR analysis methods using a common dataset. To this end, as part of the Microbial Bioinformatics Hackathon and Workshop 2021, a 'gold standard' reference genomic and simulated metagenomic dataset was generated containing raw sequence reads mapped against their corresponding reference genome from a range of 174 potentially pathogenic bacteria. These datasets and their accompanying metadata are freely available for use in benchmarking studies of bacteria and their antimicrobial resistance genes and will help improve tool development for the identification of AMR genes in complex samples.

Tatajuba: exploring the distribution of homopolymer tracts.

Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialized bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and help predict their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source licence GNU GPLv3.

SARS-CoV-2 variants of concern dominate in Lahore, Pakistan in April 2021.

The SARS-CoV-2 pandemic continues to expand globally, with case numbers rising in many areas of the world, including the Indian sub-continent. Pakistan has one of the world's largest populations, of over 200 million people and is experiencing a severe third wave of infections caused by SARS-CoV-2 that began in March 2021. In Pakistan, during the third wave until now only 12 SARS-CoV-2 genomes have been collected and among these nine are from Islamabad. This highlights the need for more genome sequencing to allow surveillance of variants in circulation. In fact, more genomes are available among travellers with a travel history from Pakistan, than from within the country itself. We thus aimed to provide a snapshot assessment of circulating lineages in Lahore and surrounding areas with a combined population of 11.1 million. Within a week of April 2021, 102 samples were sequenced. The samples were randomly collected from two hospitals with a diagnostic PCR cutoff value of less than 25 cycles. Analysis of the lineages shows that the Alpha variant of concern (first identified in the UK) dominates, accounting for 97.9 % (97/99) of cases, with the Beta variant of concern (first identified in South Africa) accounting for 2.0 % (2/99) of cases. No other lineages were observed. In depth analysis of the Alpha lineages indicated multiple separate introductions and subsequent establishment within the region. Eight samples were identical to genomes observed in Europe (seven UK, one Switzerland), indicating recent transmission. Genomes of other samples show evidence that these have evolved, indicating sustained transmission over a period of time either within Pakistan or other countries with low-density genome sequencing. Vaccines remain effective against Alpha, however, the low level of Beta against which some vaccines are less effective demonstrates the requirement for continued prospective genomic surveillance.

Genomic epidemiology and the role of international and regional travel in the SARS-CoV-2 epidemic in Zimbabwe: a retrospective study of routinely collected surveillance data.

BACKGROUND: Advances in SARS-CoV-2 sequencing have enabled identification of new variants, tracking of its evolution, and monitoring of its spread. We aimed to use whole genome sequencing to describe the molecular epidemiology of the SARS-CoV-2 outbreak and to inform the implementation of effective public health interventions for control in Zimbabwe. METHODS: We performed a retrospective study of nasopharyngeal samples collected from nine laboratories in Zimbabwe between March 20 and Oct 16, 2020. Samples were taken as a result of quarantine procedures for international arrivals or to test for infection in people who were symptomatic or close contacts of positive cases. Samples that had a cycle threshold of less than 30 in the diagnostic PCR test were processed for sequencing. We began our analysis in July, 2020 (120 days since the first case), with a follow-up in October, 2020 (at 210 days since the first case). The phylogenetic relationship of the genome sequences within Zimbabwe and global samples was established using maximum likelihood and Bayesian methods. FINDINGS: Of 92 299 nasopharyngeal samples collected during the study period, 8099 were PCR-positive and 328 were available for sequencing, with 156 passing sequence quality control. 83 (53%) of 156 were from female participants. At least 26 independent introductions of SARS-CoV-2 into Zimbabwe in the first 210 days were associated with 12 global lineages. 151 (97%) of 156 had the Asp614Gly mutation in the spike protein. Most cases, 93 (60%), were imported from outside Zimbabwe. Community transmission was reported 6 days after the onset of the outbreak. INTERPRETATION: Initial public health interventions delayed onset of SARS-CoV-2 community transmission after the introduction of the virus from international and regional migration in Zimbabwe. Global whole genome sequence data are essential to reveal major routes of spread and guide intervention strategies. FUNDING: WHO, Africa CDC, Biotechnology and Biological Sciences Research Council, Medical Research Council, National Institute for Health Research, and Genome Research Limited.

Large-scale sequencing of SARS-CoV-2 genomes from one region allows detailed epidemiology and enables local outbreak management.

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organizations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak analyses. Overall, 1565 positive samples (172 per 100 000 population) from 1376 cases were evaluated; for 140 cases between two and six samples were available providing longitudinal data. This represented 42.6 % of all positive samples identified by hospital testing in the region and encompassed those with clinical need, and health and care workers and their families. In total, 1035 cases had genome sequences of sufficient quality to provide phylogenetic lineages. These genomes belonged to 26 distinct global lineages, indicating that there were multiple separate introductions into the region. Furthermore, 100 genetically distinct UK lineages were detected demonstrating local evolution, at a rate of ~2 SNPs per month, and multiple co-occurring lineages as the pandemic progressed. Our analysis: identified a discrete sublineage associated with six care facilities; found no evidence of reinfection in longitudinal samples; ruled out a nosocomial outbreak; identified 16 lineages in key workers which were not in patients, indicating infection control measures were effective; and found the D614G spike protein mutation which is linked to increased transmissibility dominates the samples and rapidly confirmed relatedness of cases in an outbreak at a food processing facility. The large-scale genome sequencing of SARS-CoV-2-positive samples has provided valuable additional data for public health epidemiology in the Norfolk region, and will continue to help identify and untangle hidden transmission chains as the pandemic evolves.

CoronaHiT: high-throughput sequencing of SARS-CoV-2 genomes.

We present CoronaHiT, a platform and throughput flexible method for sequencing SARS-CoV-2 genomes (≤ 96 on MinION or > 96 on Illumina NextSeq) depending on changing requirements experienced during the pandemic. CoronaHiT uses transposase-based library preparation of ARTIC PCR products. Method performance was demonstrated by sequencing 2 plates containing 95 and 59 SARS-CoV-2 genomes on nanopore and Illumina platforms and comparing to the ARTIC LoCost nanopore method. Of the 154 samples sequenced using all 3 methods, ≥ 90% genome coverage was obtained for 64.3% using ARTIC LoCost, 71.4% using CoronaHiT-ONT and 76.6% using CoronaHiT-Illumina, with almost identical clustering on a maximum likelihood tree. This protocol will aid the rapid expansion of SARS-CoV-2 genome sequencing globally.

Taxonomic resolution of the ribosomal RNA operon in bacteria: implications for its use with long-read sequencing.

DNA barcoding through the use of amplified regions of the ribosomal operon, such as the 16S gene, is a routine method to gain an overview of the microbial taxonomic diversity within a sample without the need to isolate and culture the microbes present. However, bacterial cells usually have multiple copies of this ribosomal operon, and choosing the 'wrong' copy could provide a misleading species classification. While this presents less of a problem for well-characterized organisms with large sequence databases to interrogate, it is a significant challenge for lesser known organisms with unknown copy number and diversity. Using the entire length of the ribosomal operon, which encompasses the 16S, 23S, 5S and internal transcribed spacer regions, should provide greater taxonomic resolution but has not been well explored. Here, we use publicly available reference genomes and explore the theoretical boundaries when using concatenated genes and the full-length ribosomal operons, which has been made possible by the development and uptake of long-read sequencing technologies. We quantify the issues of both copy choice and operon length in a phylogenetic context to demonstrate that longer regions improve the phylogenetic signal while maintaining taxonomic accuracy.

Nanoneedle-Mediated Stimulation of Cell Mechanotransduction Machinery.

Biomaterial substrates can be engineered to present topographical signals to cells which, through interactions between the material and active components of the cell membrane, regulate key cellular processes and guide cell fate decisions. However, targeting mechanoresponsive elements that reside within the intracellular domain is a concept that has only recently emerged. Here, we show that mesoporous silicon nanoneedle arrays interact simultaneously with the cell membrane, cytoskeleton, and nucleus of primary human cells, generating distinct responses at each of these cellular compartments. Specifically, nanoneedles inhibit focal adhesion maturation at the membrane, reduce tension in the cytoskeleton, and lead to remodeling of the nuclear envelope at sites of impingement. The combined changes in actin cytoskeleton assembly, expression and segregation of the nuclear lamina, and localization of Yes-associated protein (YAP) correlate differently from what is canonically observed upon stimulation at the cell membrane, revealing that biophysical cues directed to the intracellular space can generate heretofore unobserved mechanosensory responses. These findings highlight the ability of nanoneedles to study and direct the phenotype of large cell populations simultaneously, through biophysical interactions with multiple mechanoresponsive components.

Species Tree Estimation from Genome-Wide Data with guenomu.

The history of particular genes and that of the species that carry them can be different for a variety of reasons. In particular, gene trees and species trees can differ due to well-known evolutionary processes such as gene duplication and loss, lateral gene transfer, or incomplete lineage sorting. Species tree reconstruction methods have been developed to take this incongruence into account; these can be divided grossly into supertree and supermatrix approaches. Here we introduce a new Bayesian hierarchical model that we have recently developed and implemented in the program guenomu. The new model considers multiple sources of gene tree/species tree disagreement. Guenomu takes as input posterior distributions of unrooted gene tree topologies for multiple gene families, in order to estimate the posterior distribution of rooted species tree topologies.

Infinitely long branches and an informal test of common ancestry.

BACKGROUND: The evidence for universal common ancestry (UCA) is vast and persuasive. A phylogenetic test has been proposed for quantifying its odds against independently originated sequences based on the comparison between one versus several trees. This test was successfully applied to a well-supported homologous sequence alignment, which was however criticized on the basis of simulations showing that alignments without any phylogenetic structure could mislead its conclusions. RESULTS: Here we present a simplified version of this same counterexample, which can be interpreted as a tree with arbitrarily long branches, and where the UCA test fails again. We also present another case whereby any sufficiently similar alignment will favour UCA irrespective of the true independent origins for the sequences. Finally, we present a class of frequentist tests that perform better than the purportedly formal UCA test. CONCLUSION: Despite claims to the contrary, we show that the counterexamples successfully detected a drawback of the original UCA test, of relying on sequence similarity. In light of our own simulations, we therefore conclude that the UCA test as originally proposed should not be trusted unless convergence has already been ruled out a priori.

SimPhy: Phylogenomic Simulation of Gene, Locus, and Species Trees.

We present a fast and flexible software package--SimPhy--for the simulation of multiple gene families evolving under incomplete lineage sorting, gene duplication and loss, horizontal gene transfer--all three potentially leading to species tree/gene tree discordance--and gene conversion. SimPhy implements a hierarchical phylogenetic model in which the evolution of species, locus, and gene trees is governed by global and local parameters (e.g., genome-wide, species-specific, locus-specific), that can be fixed or be sampled from a priori statistical distributions. SimPhy also incorporates comprehensive models of substitution rate variation among lineages (uncorrelated relaxed clocks) and the capability of simulating partitioned nucleotide, codon, and protein multilocus sequence alignments under a plethora of substitution models using the program INDELible. We validate SimPhy's output using theoretical expectations and other programs, and show that it scales extremely well with complex models and/or large trees, being an order of magnitude faster than the most similar program (DLCoal-Sim). In addition, we demonstrate how SimPhy can be useful to understand interactions among different evolutionary processes, conducting a simulation study to characterize the systematic overestimation of the duplication time when using standard reconciliation methods. SimPhy is available at https://github.com/adamallo/SimPhy, where users can find the source code, precompiled executables, a detailed manual and example cases.

A Bayesian Supertree Model for Genome-Wide Species Tree Reconstruction.

Current phylogenomic data sets highlight the need for species tree methods able to deal with several sources of gene tree/species tree incongruence. At the same time, we need to make most use of all available data. Most species tree methods deal with single processes of phylogenetic discordance, namely, gene duplication and loss, incomplete lineage sorting (ILS) or horizontal gene transfer. In this manuscript, we address the problem of species tree inference from multilocus, genome-wide data sets regardless of the presence of gene duplication and loss and ILS therefore without the need to identify orthologs or to use a single individual per species. We do this by extending the idea of Maximum Likelihood (ML) supertrees to a hierarchical Bayesian model where several sources of gene tree/species tree disagreement can be accounted for in a modular manner. We implemented this model in a computer program called guenomu whose inputs are posterior distributions of unrooted gene tree topologies for multiple gene families, and whose output is the posterior distribution of rooted species tree topologies. We conducted extensive simulations to evaluate the performance of our approach in comparison with other species tree approaches able to deal with more than one leaf from the same species. Our method ranked best under simulated data sets, in spite of ignoring branch lengths, and performed well on empirical data, as well as being fast enough to analyze relatively large data sets. Our Bayesian supertree method was also very successful in obtaining better estimates of gene trees, by reducing the uncertainty in their distributions. In addition, our results show that under complex simulation scenarios, gene tree parsimony is also a competitive approach once we consider its speed, in contrast to more sophisticated models.